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Abstract 

Employing the QCD light-cone sum rule approach we calculate the B — > tttt 
hadronic matrix element of the current-current operator with c quarks in the pen- 
guin topology ("charming penguin"). The dominant contribution to the sum rule 
is due to the c-quark loop at short distances and is of 0(a s ) with respect to the 
factorizable B — > tttt amplitude. The effects of soft gluons are suppressed at least 
by 0(a s m^ ). Our result indicates that sizable nonperturbative effects generated 
by charming penguins at finite mj, are absent. The same is valid for the penguin 
contractions of the current-current operators with light quarks. 
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1. Charmless two-body hadronic B decays, such as B — > tttt, are promising sources of 
information on CP-violation in the 6-flavour sector. The task of unfolding CKM phases 
from the decay observables is challenged by our limited ability to calculate the relevant 
hadronic interactions. One has to resort to approximate methods based on the expansion 
in the inverse 6-quark mass. The long-distance effects are then parametrized in terms of 
process-independent characteristics: heavy-light form factors, hadronic decay constants 
and light-cone distribution amplitudes (DA). For example, in the QCD factorization 
approach pQ the exclusive S-decay amplitudes in the m;, — > oo limit are expressed in 
terms of the factorizable part and calculable 0(a s ) nonfactorizable corrections. For 
phenomenological applications it is important to investigate the subleading effects in 
the decay amplitudes suppressed by inverse powers of m&. Especially interesting are 
"soft" nonfactorizable effects, involving low-virtuality gluons and quarks, not necessarily 
accompanied by an a s -suppression. 

Quantitative estimates of nonfactorizable contributions, including the power-suppressed 
ones can be obtained [2] using the method of QCD light-cone sum rules (LCSR) [3]. In 
particular, the 0(l/m b ) soft-gluon corrections have been calculated |2] for the B — > rnr 
matrix elements with the emission topology. Furthermore, the LCSR estimate of the 
chromomagnetic dipole operator (gluonic penguin) contribution to B — > Tin was ob- 
tained in jlj. In this case the soft-gluon contribution, being suppressed as I /ml, at finite 
rrib is of the magnitude of the 0(a s ) hard-gluon part. However, the nonfactorizable 
corrections are nonuniversal, depending on the effective operators and quark topologies 
involved in a given B decay channel, therefore these corrections have to be investigated 
one by one. 

Among the most intriguing effects in charmless B decays are the so called "charming 
penguins". The c-quark pair emitted in the b — > ccd(s) decay propagates in the environ- 
ment of the light spectator cloud and annihilates to gluons, the latter being absorbed in 
the final charmless state. In this, so called BSS-mechanism [S] the intermediate cc loop 
generates an imaginary part, contributing to the final-state strong rescattering phase. 
In QCD factorization approach [I], charming penguins are typically small, being a part 
of the 0(a s ) nonfactorizable correction to the B — ► Tin amplitude. On the other hand, 
fits of two-body charmless B decays do not exclude substantial 0(l/m&) nonperturba- 
tive effects of the charming-penguin type [Hj. It is therefore important to investigate by 
independent methods the effects generated by c-quark loops in charmless B decays. 

In this letter we report on the LCSR estimate of the penguin topology contributions 
to the B — > mi amplitude. We start with calculating the c-quark part of this effect 
(charming penguin). Later on, we extend the sum rule analysis to the penguin contrac- 
tions of the M-quark current-current operators. 

2. As a study case we choose the B° — > 7r + 7r~ channel. The decay amplitude is given 
by the hadronic matrix element (7r + 7r~|i/ c fj|I? ) of the effective weak Hamiltonian [7] 

H « = % { E W + c 2 oi) - v tb v; d (j2 c i°i + c%A* ) } , (i) 

v ^ p=u,c \ i=3 / ) 

where 0\ = (dT ^(pT^b) and O v 2 = (pT '^(dT^b) are the current-current operators 
(p = u,c and T M = 7^(1 — 75)), O3-10 are the penguin operators, and 0$ g is the 
chromomagnetic dipole operator. Each operator entering Eq. (0) contributes to the 
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Figure 1: Diagrams corresponding to the 0(a s ) penguin contractions in the correlation 
function (HJ). Only the diagrams contributing to the sum rule are shown. The square 
denotes the four-quark operator (D\. 

B — > mi decay amplitude with a number of different contractions of the quark lines 
(topologies) f5|. In what follows we will mainly concentrate on the operator 0\. For 
convenience we decompose this operator 0\ = \0^ + 202, extracting the colour-octet 
part O c 2 = (cr^f c)(dT^ b) . 

To derive LCSR for the B — > ixix hadronic matrix element of 0\ we follow the pro- 
cedure which is described in detail in j2J E] . One starts from introducing the correlation 
function: 

Ff S) = i 2 Jd A xe-^>J d%e^y(0\T{j ( -\y)dmji B \x)}\n-(q)) 

= (p-k) a F( Sl ,s 2 ,P 2 ) + ..., (2) 

where jj^ = u'-faj^d and = imjry^d are the quark currents interpolating pion and B 
meson, respectively. In the above, ellipses denote the Lorentz-structures which are not 
used. Only the colour-octet part of 0\ needs to be taken into account. The contributions 
of 0\ contain at least a two-gluon annihilation of the color-neutral c-quark pair and have 
to be considered within higher-order corrections. The correlation function (J2J) depends 
on the artificial four-momentum k allowing one to avoid overlaps of the 6-flavoured and 
light-quark states in the dispersion relation. We also choose p 2 = k 2 = for simplicity 
and adopt the chiral limit q 2 = m 2 = 0. Thus, the invariant amplitude F is a function 
of three variables 8\ = (p — k) 2 , s 2 — (p — o) 2 and P 2 = (p — q — k) 2 . 

The next step is to calculate F(si, s%, P 2 ) at large spacelike Si, s%, P 2 employing the 
operator product expansion (OPE) near the light-cone. The corresponding diagrams of 
0(a s ) are shown in Fig. ^ They contain a c-quark loop, which involves a well known 
function of m 2 . The divergence of the quark loop is absorbed in the renormalization 
of the QCD penguin operators C?3_6- The finite contribution of the loop is formally of 
the next-to-leading order and depends on the renormalization scale \x. This dependence 
is compensated by the Wilson coefficients of the QCD penguin operators so that the 
remaining \i dependence is weak, since it is of 0(a 2 ln /i). Furthermore, the finite piece 
of the loop contains scheme- dependent constants; we use the NDR scheme jTj. 

The sum rule is derived from the dispersion relations for F in the variables s± and 
S2- Therefore the calculation of the two- loop diagrams in Fig. can be considerably 
simplified if one starts from their imaginary parts obtained by employing the Cutkosky 
rule and replacing the propagators by 5-functions. For the diagram in Fig. it is 
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sufficient to calculate its imaginary part in s\. Furthermore, according to the procedure 
explained in [2] we have to consider the dispersion relation in Si in the quark-hadron 
duality interval, < si < Sq <C m\ and simultaneously at |P 2 | ~ m 2 B . Hence, one can 
safely neglect all 0(si/P 2 ) contributions in Im sl F(si, s 2 , P 2 ). In this approximation we 
obtain for the diagram in Fig. ^i: 



Im Sl Ff 9d (s 1 ,s 2 ,P 2 ) 



a s C F fnml 
8tt 2 



du 



mi - us 2 
o o 



dz I(zuP 2 , m c ) 



2rrib 

S 2 \ / (p a {u) U(p' a (u) 



xP 2 ^z(l-z)<p„(u) + (l-z) 
+ (2z 



P 2 



6 



2z + Jl ) u M u ) 



(3) 



In the above, <p n and <p p , ip a are the pion DA of twist 2 and 3, respectively, the latter 
are normalized by = m 2 /(m u + m^) nonvanishing in the chiral limit; we use the same 
standard definitions as in [I]. Finally, f' a (u) = dip a (u)/du and !(...) is the c-quark loop 
function in the NDR scheme: 



dxxil — x) In 



x(l — x)f 



(4) 



Similarly, for the diagram in Fig. QJd it is sufficient to calculate its imaginary part in s 2 . 
The result reads: 



Im S2 F^\ Sl ,s 2 ,P 2 



a s C F f^m 2 h fs 2 -m 



16tt 2 



2\ 2 



S-2 



du 



uP 2 + US\ 
o o 



dzl c (—uz(s 2 — ml)) 



x < (P 2 - Sl - s 2 )z <p r (u) + i^-y 2 - si + 3s 2 )v P (u) 



12m& 



2 I ^ 31 82 + 2 _ P 1 . 31 (P 2 -Si- 82) I r -„( « ) 



uP 2 + itSi 



+ (3P 2 - 3 Sl + s 2 V» 



(5) 



The remaining diagrams not shown in Fig. Q with gluons attached to the virtual b and 
d lines, do not contribute to the sum rule because their double imaginary parts vanish 
inside the duality regions 0< s± < Sq, m\ < s 2 < Sq. 

Eqs. (JHJ) and (JHJ provide the leading, 0(a s ) and twist 2 and 3 answer for the cor- 
relation function. The contributions of the higher-twist pion DA's to F are suppressed 
by inverse powers of large squares of external momenta, reflecting the power counting 
in the light-cone OPE. Accordingly, in this paper, we neglect small effects of the twist-4 
quark- ant iquark DA's in the diagrams of Fig. 

3. So far we have taken into account the hard-gluon emission off the c-quark loop. 
The most intriguing and difficult problem concerning charming penguins is the effect 
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Figure 2: Diagrams containing the quark- antiquark-gluon DA's 

of the soft (low-virtuality) gluons coupled to the c-quark loop. Within the sum rule 
approach this problem has to be addressed in the context of the light-cone OPE of the 
correlation function. Together with u and d quarks, the on-shell gluons emitted at short 
distances form multiparticle DA's of the pion. Importantly, the contributions of these 
DA's to the sum rule, being of higher twist, are suppressed OH] by inverse powers of the 
heavy mass scale with respect to the contributions of 2-particle quark-antiquark DA's of 
lower twists. Therefore, for a sum rule estimate of the soft-gluon effects it is sufficient to 
consider diagrams with one "constituent" gluon i.e., diagrams involving quark- antiquark- 
gluon DA's of the pion. For the gluonic penguin operator, this soft-gluon effect turns 
out to be important [4], because in the quark-antiquark-gluon term of the sum rule the 
I /ml suppression is compensated by the absence of a s . In the correlation function (J2J) 
the situation is different. The diagram with one gluon shown in FigEt vanishes due to 
the current conservation in the c-quark loop. Nonvanishing terms with the three-particle 
DA's emerge from the diagrams, containing at least one hard gluon in addition to the 
on-shell gluon. One of these diagrams is shown in Fig|2b- Their complete calculation is 
a difficult task. From the studies of the b — > s'j matrix elements of 0i,2, where similar 
diagrams with an on-shell photon and virtual gluon have been calculated |Hj , we conclude 
that the contribution of Fig. Eb and similar diagrams not only contain a s but also have 
an additional 0(l/m|) suppression with respect to the diagrams in Fig. ^ 

The next nonvanishing contribution to the expansion of the coloured c-quark loop 
near the light-cone contains a derivative of the gluon field D u G a ^ u which can be further 
reduced to the light-quark pair due to QCD equation of motion. The whole effect has 
a short- distance nature, similar to the formation of the QCD penguin operators. The 
resulting diagram shown in Fig. has to be calculated in terms of the pion four-quark 
DA's and is beyond the approximation adopted here and in [4]. Simple dimension count- 
ing yields for this diagram a suppression factor of 0(l/m|). The c-loop factor in this case 
reduces to \n.{m 2 c / fi 2 ) plus a scheme- dependent constant. The /i- and scheme-dependence 
are compensated by the contributions of QCD penguin operators, since the c-quark loop 
momenta larger than fi are included in the short- distance coefficients of these operators. 

The diagram with two gluons emitted from the c-quark loop (Fig. Eb) is also not 
included in our calculation because it contains DA's with multiplicity larger than three. 
Here both color- neutral and color-octet parts of 0\ (that is, both 2 and 2 ) con- 
tribute. For simplicity we concentrate on 2 where only the axial-vector part of the 
c-quark current is relevant. Since all light degrees of freedom in this diagram en- 
ter DA of the massless pion (q 2 = 0), it is clear that the four- momentum flowing 
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(a) (b) (c) 

Figure 3: (a),(b) Diagrams corresponding to the multiparticle DA's of the pion and 
arising from the expansion of the c-quark loop in the external gluon field; (c) the four- 
quark "soft" contribution in the case of the operators 0" 2 . 



through the c-quark current has to be lightlike, i.e. the c-loop remains far off-shell 
(with a typical virtuality scale of 2m c independent of mj. Hence, one can use the 
local expansion of the c-quark axial- vector current in the gluon field obtained in [TUj : 



COW -> g 2 J(1927r 2 m 2 c ) 



The suppression factor with 



d^{G a pu G apv ) - A(D a G a ' ua )G a ^^ 

respect to the diagrams in Fig. ^obtained by counting dimensions of the resulting oper- 
ators (DA's) is for this contribution at least 0(1/ (m 2 m 2 )). The presence of ln(m c ) and 
m~ 2 in the contributions of Fig|Hk and 3b respectively, indicates that at m c — > these 
terms are divergent. In other words, if c quarks are replaced with the light quarks, e.g., 
in the case of 0", the light-quark pair propagates at long distances, that is, belongs to 
the pion four-quark DA. The corresponding diagram is the one shown in Fig. 

As already mentioned, we neglect the contributions of four-quark DA's stemming 
from the matrix elements of the type (0 | u{xi)q(x 2 )q(x 3 )d(x i )\'w) (xi on the light-cone). 
On the other hand, following jlj we take into account the factorizable parts of the 4-quark 
vacuum-pion matrix elements, extracting the configurations where one quark-antiquark 
pair forms the quark vacuum condensate, whereas the other one hadronizes into a twist 
2 and 3 pion DA. Such contributions are enhanced by the large parameter fi n . In the 
approximation adopted here, only two diagrams shown in Fig. 0] contribute to the sum 
rule (that is, have a nonvanishing contribution to the double dispersion relation in the 
duality region). Note that the quark-condensate diagram in Fig. ^jp originates from 
the 4-quark diagram in Fig. E^. The diagrams in Fig. 0] have only one loop and their 
calculation is relatively simple yielding the following result: 



+3/j, 7r m b ip p (u)-{ 




a s C F f n m b (qq) 
12tt 



P 2 - Sl 
uP 2 + USi 



du 



[mi - us 2 )s 1 



I(uP +usi,m c ) 



US 2 



S-2 



m 2 — us 2 J ml — us 2 



<p<t( u )-<pM 



P 2 - 3^2 - Sl ( P 2 - Sl , 

</M«) - 3y? p (w) H Lp a [u 

2 V Qs 2 




■ (6) 



Multiparticle contributions which are factorized in the condensates of higher dimension 
are not taken into account. We also neglect the quark-condensate contributions of the 
type (qq)(0 \ u(xi)G^ u (x 2 )d(x3) \ ir) arising from the diagram in Fig. Ob after applying 
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(a) (b) 

Figure 4: Diagram corresponding to the factorizable J^-quark contribution to the correla- 
tion junction. 



the QCD equation of motion to the derivatives of G" . These terms are suppressed at 
least by 0{l/m 2 ) with respect to the diagrams in Fig. |3J 

Summarizing, we do not find significant contributions involving soft gluons in the 
OPE of the correlation function (J2J). The dominant effect arises from the c-quark loop 
annihilation into hard gluons (Fig. QJ. In addition, there is the quark-condensate contri- 
bution (Fig. EJ) which we consider a natural upper limit for all neglected contributions of 
multiparticle DA's. Importantly, the presence of the new intermediate scale m c <C in 
the correlation function does not noticeably enhance the charming penguin contribution: 
in the leading terms of OPE m c enters logarithmically, and the inverse powers of m c 
appear only in the subleading suppressed terms. 

4. Having calculated the relevant contributions to the correlation function we are now 
in a position to obtain the sum rule for the hadronic matrix element A^ *' (B® — > = 
(n^ (p)n + (—q)\C>2\B^(p — q)) . The derivation of LCSR explained in detail in E] 
includes: use of the dispersion relations for the invariant amplitude F(s\, s 2 , P 2 ) in the 
variables si and s 2 ; employing quark-hadron duality in both pion and B meson channels 
with the threshold parameters Sq and Sq, respectively; the Borel transformations and 
the continuation from large spacelike P 2 to large timelike P 2 = m 2 B . As a result, LCSR 
is given by the following expression: 

Uf B A^\B° d ^^-)e^m 
= Jds 2 e~^ M2 4j d Sl e-^ M hm S2 lm Sl F( Sl ,s 2 ,P 2 ) 

where M\ and M 2 are the Borel parameters in the pion and 5-meson channels, re- 
spectively. Using Eqs. ©,© and (0), we obtain from Eq. (JZJ) the "charming penguin" 
hadronic matrix element at 0(a s ) and with twist 2 and 3 accuracy : 

A^\B° d - 7T+7T-) = A FigA + A m , (8) 




7 



where 



AFiq.l 




dse S ' M * \ m\U fdu (m |_i^ )/M | 



bTb J u 

B 




+(l — z)—^- iup p (u)[2z + 



um 



B 



/dz I(—zmlulu) — 
u 



Ml) 1 + 



3 6 
3mg\ , ^(1) 



2z - 



B 



+ 



5m| 



(9) 



is the contribution of the diagrams in Fig. ^ and 

i , 
'Gir b { Um h ) Wlfs) J u* 6 P 



A 



um B ,m 2 c ) 



2ip 7T (u) + — 3u(f p (u) + 

m b \ 3 



+/(0,m* 



3 — H 3uv? p (m) + ^(m) 



m 6 



6m| 



(10) 



is the factorizable 4-quark contribution. In the above, Uq = m 2 /sQ. The sum rule (|1U|) is 
obtained at finite but we neglect numerically very small corrections of order s^/m^. 
Note that the c-quark loop factor I(zm B u, m 2 ) in the sum rule originating from the 
Fig. diagram has a complex phase which is generated by the analytic continuation 
in P 2 and has to be associated with the strong-interaction phase (in the quark-hadron 
duality approximation). The BSS mechanism ^5, is thus recovered in LCSR. 

To compare our result with QCD factorization we investigate the heavy-quark mass 
limit of the sum rule (jHJ), expanding all heavy-mass dependent quantities in powers of 
mt,. At m>b — > oo only the diagram of Fig. contributes. In this limit duality interval of 
the M-integration in Eq. © reduces to the end-point region allowing one to put u — > 1 in 
the c-quark loop factor. After that, the integral over u multiplying this factor reproduces 
LCSR for the B — > n form factor ^T]. To see that one has to use the relation |12J following 
from the QCD equation of motion : mp p {u) + <p a (u)/3 — wp' a {u)/§ = (neglecting the 
quark- ant iquark gluon DA's). After these transformations we obtain : 



7T 7T 



71 



dzz(l-z)I(zmt,ml) ) ^ff VZ) fL(0) [LCSR) ) ■ 



,{LCSR)\ 



mi,— >oo 



In the above, we have denoted the integral over s in Eq. (jHJ) by fi svz ^ indicating that it 
coincides with the leading-order SVZ sum rule for /„. |13j . Furthermore we have replaced 
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the LCSR B — > tt form factor by an m^-independent effective form factor /£ vr (0)( iC ' 5 ' R ) = 

fBn(®)^ LCSR ^ / m b^ 2 ■ Thus, the charming penguin matrix element at mj — > oo factorizes 
into the c-loop integral and the factorizable B — > 7r7r amplitude equal to the matrix 
element of 0\ in the emission topology: A^\B d — ► Tr" 1 "^"^-^ = iy/mbf n f^(Q). 
Furthermore, since the asymptotic twist-2 pion DA has the form ^(z) = 6^(1 — z), the 
limiting expression (jllj) reproduces the c-quark penguin-loop contribution to the B — ► 7T7t 
amplitude in the QCD factorization pQ at the twist 2 level. For brevity, we omit a more 
detailed discussion of various l/m& terms which at the end lead to a numerical deviation 
of the LCSR (JBJ) from the — > oo limit. We only note that the contribution of the 
diagram in Fig. which does not have a counterpart in QCD factorization is suppressed 
as 1/m 2 . Strictly speaking this diagram turns out to be beyond the adopted accuracy 
since in LCSR we have neglected other small 0(1 /ml) terms. 

We conclude this section with a comment on the penguin contractions of the current- 
current operator 0\. The corresponding sum rule for the hadronic matrix element 

Af^(B° d -»■ 7T+7T-) = (n-n + \0%\B%) P is easily obtained from LCSR © by putting 
m c — ► everywhere except in the second term of the quark condensate contribution 
(jlOp . As we already discussed this term emerges from the factorization of the diagram in 
Fig. Ok and in the case of the light quarks it has to be absorbed in the four-quark DA's 
(Fig. Ofc). This four-quark DA contribution does not contain any a s suppression. How- 
ever, counting the dimensions we find that it is power-suppressed, at least by 0(l/m\). 

5. For a numerical estimate of the charming penguin in B — > tttt decay we calcu- 
late the ratio of the sum rule (JBJ) to the factorizable amplitude A^ X \B^ — > ix + n~) = 

ro*>,™ + n AWXm -»• 7T+7T-) 2A^(B° d ^7V + 7r-) , x 

We use the same input as in [4- and in the numerical analysis of the LCSR for the 
B — » 7r form factor (TJ. The following parameters are taken for the pion channel: 
f w = 132 MeV, = 0.7 GeV 2 , M\ = 0.5-1.5 GeV 2 ; and for the B channel: m b = 
4.7 ± 0.1 GeV (the one-loop pole mass), s% = 35 =f 2 GeV 2 , M | = 8-12 GeV 2 . The 
normalization scale adopted for the pion DA's and a s is //& = a/ m 2 ^ — m 2 ~ 2.4 GeV. 
For the c quark mass we take m c = 1.3 ±0.1 GeV, and for the quark condensate density 
(qq)(lGeV) = -(240 ±10 MeV) 3 , or, equivalently, ^(lGeV) = 1.59 ±0.2 GeV. All pion 
light-cone DA's are taken in the asymptotic form: ip n (u) = (f a (u) = 6u(l—u), f p (u) = 1. 
To get a feeling how the nonasymptotic form of the pion DA influences our result we 
also recalculate Eq. (fT2j) using a model for tp^u) with a nonzero second Gegenbauer 
coefficient a 2 (lGeV) = 0.4. With the above input and adding the uncertainties caused 
by the variation of all parameters linearly, we get the following range: 



r 



(Of) 



(B° d -> tt+tO = [-(0.29 -T- 0.56) - (1.3 ± 1.6)z] • 10" 2 . (13) 



The corresponding estimate for the penguin contraction of the operator 0\ is: 

rf f} (B° d -> tt+tt") = [(0.09 -T- 0.21) - (1.6 ± 2.1)i] ■ 10" 2 . (14) 

The penguin-topology contributions turn out to be very small, not larger than the 
other nonfactorizable corrections in B — > tttt, however, as we shall see in the next section, 



9 



they add up to a noticeable effect in the CP asymmetry. The contribution of the quark 
condensate term in the sum rule to the real parts of Eqs. (j!3|) and (fT^j) is relatively small, 
but constitutes about 50% of the imaginary parts. 

6. It is well known that the penguin contributions play a key role in the direct and 
mixing-induced CP violation in B — > tttt. The time-dependent CP-asymmetry is given 
by 

Acp{B% -> 7i + 7i-)(t) = a^cos(AM d t) + a^sm{AM d t) , (15) 

where ag£ = (1 - |£| 2 )/(1 + |£| 2 ) and a^ x = (2Im£)/(l + |£| 2 ), with f = e - 2i{p+ ^(l + 
Re^)/{1 + Re~^) and R = -P/(R b T). Here T is the contribution to the B -> tttt 
amplitude proportional to V u bV* d = \V u bV* d \e~ n and contains also the penguin contrac- 
tion of the current- current operator 0\. The remaining amplitude P is proportional to 
VcbV* d . The factor Rb = \ V u b\ \V ud \/(\ V c b\ \V cd \) as usual defines one side of the unitarity 
triangle. We take R b = 0.39 ± 0.04. 

Strong phases originate from both T and P; thus we have T = |T|e J<5T and P = \P\e lSp 
and 





-2\R\ 


sin(<5p — St) sin 7 


1 - 2\R\ 


cos(5p — 5t) cos 7 + 


R\ 


2 



(16) 



In the previous sections we have calculated the contributions arising from the penguin 
contractions of the operators 0\ and 0\. In addition to these we have also to take 
into account the tree and penguin contributions of the penguin operators £?3_6 and the 
gluonic penguin contribution of the dipole operator (D% g . For the latter we use the LCSR 
result from |1J. The electroweak penguin contributions to B — > tttt are color-suppressed 
and very likely negligible. 

Apart from the penguin contractions one also has to include the hard nonfactorizable 
0(a s ) and soft l/m& corrections (in the emission topology) to both T and P, obtained, 
respectively in the QCD factorization approach pQ and from LCSR 0. It turns out 
that their contributions to the CP-asymmetry are negligibly small. To this end, we can 
numerically study a^?p using the penguin contributions calculated from LCSR at finite m;, 
and compare the result to the infinite-mass limit that agrees with the QCD factorization 
prediction PQ. The result for function of 7 is shown in Fig. In this calculation 

the Wilson coefficients in H e ff are taken at the same scale //& ~ mb/2 as used in LCSR. 
Let us emphasize that this is not yet the final prediction because the annihilation effects 
in B — ► tttt are not included in this analysis. The main aim of this numerical illustration 
is to demonstrate that there is a difference between the finite m>b and the m& — > 00 result 
and that this difference could be sizable (within the uncertainty of the LCSR prediction). 

7. Concluding, we estimated the B — » tttt hadronic matrix elements with penguin 
topology for the current- current operators with c and u quarks using the LCSR approach. 
The main contribution to the sum rule stems from the 0(a s ) quark loop annihilating 
to a hard gluon. This contribution determines the strong phase of the hadronic matrix 
element justifying the use of the (perturbative) BSS mechanism. The soft-gluon effects, 
which in the sum rule approach correspond to multiparticle pion DA's, are suppressed, 
at least by 0{a s /ra b v ). In m& — > 00 limit our result agrees with the QCD factorization 
prediction for the penguin contractions. In both approaches a small value of the direct 
CP asymmetry in B^ — > tt + tt~ is expected, due to the fact that the strong phase is 
generated perturbatively. At finite m;, we predict an even smaller numerical value for this 
effect, indicating that 0(a s /m b ) corrections are important. Our result does not support 
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Figure 5: Direct CP asymmetry in B — ► titt as a function of the CKM angle 7. The 
upper curve is the result obtained for — > 00. The dark region is the LCSR result, with 
all uncertainties from the method included (uncertainties in the CKM matrix elements 
are not taken into account). The light region shows the deviation from the nib 00 limit 
result. 



models employing hadronic dispersion relations with intermediate charmed meson pairs, 
predicting sizable charming-penguin effects The method and the results of this 

paper can be used for analysing penguin effects in other charmless two-body B decays. 
The remaining problem for these decays is the estimate of the annihilation effects which 
is however a very challenging and technically difficult task even for the LCSR method. 
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